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ABSTRACT 

We present non-adiabatic hydrodynamic simulations of a supersonic light jet propa- 
gating into a fully ionized medium of uniform density on a scale representative of the 
narrow line region (NLR) in Seyfert galaxies with associated radio jets. In this regime 
the cooling distance of the swept up gas in the bowshock of the jet is of the same order 
as the transverse extent of the jet bowshock, as opposed to the more extreme regimes 
found for more powerful adiabatic large scale jets or the slow galactic jets which have 
been simulated previously. We calculate the emissivity for the Ha line and radio syn- 
chrotron emission. We find that the structure of the line emitting cold envelope of the 
jet cocoon is strongly dependent on the non-stationary dynamics of the jet head as 
it propagates through the ambient medium. We observe the formation of cloud-like 
high density regions which we associate with NLR clouds and filaments. We find that 
some of these clouds might be partially neutral and represent sites of jet induced star 
formation. The calculated Ha flux and the spectral line width are consistent with 
NLR observations. The simulation of the radio-optical emission with radiative cooling 
confirms the basic result of the geometric bowshock model developed by Taylor et 
al. (1989) that the start of noticeable optical line emission can be significantly offset 
from the hotspot of the radio emission. However, the time-dependent nature of the jet 
dynamics implies significant differences from their geometric bowshock model. 

Key words: galaxies: active - galaxies: jets - galaxies: kinematics and dynamics - 
galaxies: Seyfert - galaxies: individual: Mknl066 - hydrodynamics; numerical 



1 INTRODUCTION 

The study of the Narrow Line Region (NLR) in Seyfert 
galaxies is contributing substantially to the understanding 
of the central engine of active galaxies. Of particular interest 
is the apparent close association between radio continuum 
and optical line emission in this region, which has a typi- 
cal extent of 0.1-1 kpc. The link between optical and radio 
emission in the NLR manifests itself most prominently in 
a positional association between linear radio structures and 
optical line emission as well as in a correlation between line 
width and radio power (Wilson and Willis 1980, Whittle 
1985). It is widely accepted that the radio emission is re- 
lated to the outflow of relativistic plasma from the core of 
the galaxy. As is found in higher power radio galaxies and 
quasars, this outflow appears to be highly coUimated into 
jets. The radio-optical association suggests that the interac- 
tion of the jets with the interstellar medium (ISM) strongly 



influences the dynamics of the ionized gas in the NLR. It 
can be hoped that the study of this interaction will reveal 
important information on the physical conditions of the gas 
in the NLR and in the plasma jets. 

The expansion of the cocoons of large scale extragalac- 
tic jets in radio galaxies has been examined analytically by 
Scheuer (1974) and Begelman and Cioffl (1989) and numer- 
ically by Cioffl and Blondin (1992). Blandford and Konigl 
(1979) flrst suggested that the interaction of jets with BLR 
clouds could cause the radio emission and the associated 
kinematics of the BLR. The observation of symmetric dou- 
ble or triple radio sources associated with emission line gas 
in the NLR of a number of Seyfert galaxies (in particular in 
NGC 5929, Haniff, Wilson and Ward 1988, Su et al. 1996, 
and NGC 1068, Wilson and Ulvestad 1987, Gallimore et 
al. 1996) led to the development of more detailed models 
invoking the interaction of plasmons or jets with the local 
ISM. Pedlar et al. (1985) proposed an expanding plasmon 
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model for the NLR. They suggested that a high pressure 
relativistic plasma bubble expands and drives a shock wave 
sweeping up the surrounding ISM. 

Whittle et al. (1986) found a spatial separation of 
0.95 ± 0.11" between the peaks of the [Olll] 5007-A com- 
ponents in high dispersion spectroscopic observations. This 
separation is significantly smaller than the measured dis- 
tance between the radio components. Based on these obser- 
vations, Pedlar et al. (1987) proposed a variant of the plas- 
mon model similar to the radiative bowshock model sug- 
gested earlier by Wilson & Ulvestad (1987) for the radio 
structure observed in NGC 1068. This model explains the 
difference in the spatial separation between the peaks of the 
radio and the optical components. It invokes a supersonic 
motion of a jet or plasmon inducing a bowshock in the am- 
bient medium. In the frame of the head of the bowshock, 
the shocked ISM gas then moves along its surface. All the 
basic properties of an expanding plasmon model are still ap- 
plicable, but, in addition, the cooling time of the gas behind 
the bowshock allows the gas to flow a significant distance 
from the vertex before cooling to a temperature sufficiently 
low (~ lO** K) that optical lines are emitted. Since the radio 
emission is expected to peak at the vertex, there will be a 
difference in the peak positions of the optical and the ra- 
dio emission. This model was further developed by Taylor 
et al. (1989, 1992) and Ferruit et al. (1996). 

Recent Hubble Space Telescope (e.g. Capetti et 
al. 1995a,b; Bower et al. 1994, 1995) and MERLIN observa- 
tions (Pedlar et al. 1993; Kukula et al. 1996) revealed over- 
whelming new structure in the NLR of several Seyfert galax- 
ies with radio jets. Discrete clouds and diffuse line emission 
are found to be closely aligned with highly structured and 
sometimes distorted radio jets. These observations clearly 
show that the idea of a radio jet interacting with the NLR 
can account for basic features of these jets. However, the 
phenomena are far more complicated than can be accounted 
for by the simple geometric bowshock model. The full non- 
stationary character of the interaction between a jet and 
its environment has to be taken into account to explain the 
newly observed features. 

In this paper, we therefore apply high resolution ax- 
isymmetric hydrodynamic simulations to explore the time 
dependent interaction of Seyfert jets with the NLR. We 
demonstrate that the association between the radio and op- 
tical emission can be explained as a natural consequence 
of the expansion of a hot jet cocoon into the interstellar 
medium creating an envelope of dense cool gas and discrete 
emission line knots which can be associated with the narrow- 
line clouds themselves. We compare our model with obser- 
vations of the NLR Ha emission in MKN 1066, which has a 
similar linear extension (Bower et al. 1995). 

In Section ^ we estimate the properties of our cocoon 
shock model analytically. The numerical procedure of the 
simulations is briefly described in Section |^. Section ^ con- 
tains the results and their discussion. We summarize our 
conclusions in Section |^. 



2 ANALYTICAL ESTIMATES 

We consider the dense optically emitting envelope of swept- 
up ISM gas around the cocoon of a supersonic jet propa- 



gating into a uniform, ionized medium. The ionization is 
assumed to be due to photoionization from a UV source 
at the active nucleus of the galaxy. For a UV photon rate 
S ~ 10^^ s"^ similar to the one in NGC5929 (Bower et 
al. 1995) and an ambient density of lcm~^ the radius of 
the Stromgren sphere is 700 pc, which is well beyond the 
distances we consider in this paper. For an estimate of the 
initial conditions in a simulation which can be compared 
with typical observations, we assume that the envelope is a 
thin cylindrical shell of radius ~ 20 pc, length Ic = 150 pc, 
Hq luminosity Lhq ~ lO^'^ — lO'^^ergs"^, and a thickness 
dc ^ Tc. Within our model the spectral line width in the 
NLR will be similar to twice the cocoon expansion speed, 
if observed from a direction perpendicular to the jet axis. 
We therefore have typical shock velocities of 100-200 km 
s"^(e.g. MKN 1066, Bower et al. 1995). In the following, we 
estimate some characteristic quantities, such as the density 
no of the ambient ISM, the density Tie and thickness dc of 
the dense cool envelope, using this range of observed sizes 
and luminosities. 



2.1 Non-magnetic case 

We first assume that there is no magnetic field in the en- 
vironment which could limit the compression at the shock, 
which is driven perpendicularly to the jet axis into the ISM 
by the overpressured cocoon. We further assume that the 
shock is isothermal, at the temperature To = IQ^K of the 
pre-shock gas. This assumption is justifled by estimating 
the cooling distance 5 behind a shock of speed Vs ~ 150 km 
s~^ in typical environmental gas with a hydrogen density of 
no =1 cm"'', using the adiabatic compression factor of 4. 
We then have (Taylor, Dyson and Axon 1992) 



S = 



{l-c) 



4mpnoAo 

where we use a power-law cooling function of the form 



A(re) = AoT," 



(1) 



(2) 



with Ao = 4.6 • lO'^^^ergs-^cm^* and a = -0.76 at T > 
1.5T0^K. The temperature Ts immediately behind the shock 
is given by 



3_m 2 
16 ' 



(3) 



where m is the mean molecular weight and k is the Boltz- 
mann constant (we assume m — O.bmp, m-p is the proton 
mass). The cooling distance perpendicular to the shock is 
then 5 = 1.8x10^^ cm or 0.06 pc with corresponding post- 
shock temperature of Ts = 2.5x10^ K. For Us = 100 (200) 
km s~^ the cooling distance will be 2.2 (3.7) times smaller 
(larger) than for Ws=150 km s~^. All these values are much 
smaller than the radius ~ 20pc of the cocoon, justifying 
the assumption of an isothermal shock for the estimates. 

We now calculate the density ric and thickness dc of the 
compressed postshock layer of gas and the density of the 
environmental gas no. We start with the compression ratio 
of the isothermal shock, which is 



Tic 

no 



mv, 

Tfr 



1) mp VlO^Ky VlOOkms-V 



(4) 
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In terms of the volumes of the cylindrical cocoon Vz and the 
cylindrical shell of swept up gas around the cocoon K we 
also have 



/ 



2TvrcdJc 



2dc 



(5) 



Combining equations (^) and (j^ we find the thickness of 
the layer of swept up gas to be 



rc kTp 
2 ffivi 

0.083 pc- 



20pc 



(iO-ik) ( 



100 kms-i 



(6) 



Assuming pure hydrogen gas of uniform density, the Ha lu- 
minosity of this layer will be 



= 3.7x10'"' ergs 



-^f Vf ] (7) 

V 100cm- V VlOS^cmSy ^ 



where e — 3.8x IQ-^^er gcm^ s"^ is the Ha emission coeffi- 
cient at To (see Section |3.4.1 ). 

The total mass swept up from the environment into the 
shell is 



8300Mf 



(lOOcm-0 (lOSScmO 



We assume full ionization of the gas in the NLR. A column 
of swept up hydrogen gas of density ric at a distance Ic from 
the ionizing photon source with a photon flux S will stay 
fully ionized up to a thickness 



ds — l.lpc ( — — - 

^ VlOOcm-'^ 



^ lOOpc 



(l0^2s-i) 



(9) 



assuming that no photons are intercepted on their way to 
the cloud. 

In Section |^ we compare these analytic estimates with 
the results from our simulation. 



2.2 Magnetic case 

We now consider the case where the undisturbed environ- 
ment is threaded by a random magnetic field, assuming the 
components in the plane of the shock will be compressed 
with the fleld lines frozen into the gas. The effect of the 
magnetic fleld will be to limit the compression of the shocked 
gas when the magnetic pressure equals the thermal pressure. 
The compression factor fm is then given by (e.g. Dopita and 
Sutherland 1995): 



m — 



64- 



D -1/2 

ip UOOkms-iy 



Bun, 



-1/2 



(10) 



The thickness of the shocked gas layer is then 



Bun, 



-1/2 



2vc (87rm) 
0.15 pc ffinip 



3/2 



V 20pc J VlOOkms 



B«n, 



-1/2 



/iG cm-3/2 



(11) 



Values of the magnetic parameter (_B||n(,j^^^)/(^Gcm ^^^) 
are of order 1, in the range 1-10 (Dopita & Sutherland, 
1995). 

Due to the additional magnetic pressure, the thickness 
of the cold envelope will be larger and the density lower than 
in the non-magnetic case. Since the optical line emission 
depends on the square of the density, the emission received 
from a strongly magnetized envelope will be considerably 
lower than that coming from a non-magnetic envelope. 



3 NUMERICAL SIMULATION 
3.1 Initial conditions 

The determination of jet parameters like plasma speed and 
density is a classical and basically unsolved problem in ex- 
tragalactic astrophysics. For the purpose of simulating the 
the expanding envelope, the exact values of the jet speed and 
density input parameters are not important. The expansion 
of the cocoon shock is principally determined by the co- 
coon pressure, and a given pressure can be obtained from 
a range of jet parameters. For practical purposes this range 
of parameters is also somewhat restricted by the available 
computing resources. We therefore performed a small series 
of test runs varying the particle density in the ISM around 1 
cm""^, jet velocities around 10* km s~^(Bicknell et al. 1990) 
and Mach numbers around 10. A Mach number higher than 
5 is required for a jet with a noticeable cocoon (Norman, 
Smarr and Winkler 1985). From these tests we chose a set 
of parameters which compares well with the observations of 
Seyfert galaxies given in Section ^ (in particular MKN 1066) 
and which illustrates best the qualitative features of the for- 
mation and structure of the dense cocoon envelope. A more 
detailed account on the influence of varying the parame- 
ters will be given in a forthcoming paper (Steffen et al. , in 
preparation). The chosen parameters are: 



no 


= 1.5 cm"^ 


To 


= lO^K 


nj 


= 0.18 cm"^ 


Mj 


= 8.3 


'"j 


= 6.8x10* cms 


rj 


= 5.3xl0^**cm 



This yields a jet mechanical luminosity of Lj 



2 X 



lO^^ergs-^ 



3.2 Scaling and similarity 

In the absence of cooling, the jet will be evolve in a self- 
similar fashion 



lc = A' ^ 



3\ 1/5 



mno 



with ^ ~ 1, once Ic 2> {nj /n^Y^'^r 
1991). 



(12) 



2 X 10^* cm (FaUe 
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Cooling will begin to have an effect when the advance 
speed of the outer shock at the base of the jet is roughly 
lOOkms"^ (Dyson 1984), that is roughly when Ic — icooi = 
6 X 10^® cm, if the aspect ratio of the jet is 7?. = Ic/rc ~ 5. 
Thus, for our typical parameters, the jet will be self similar 
before cooling becomes important (neglecting the effects of 
turbulence in the jet). 

Once cooling becomes non-negligible, self-similarity no 
longer holds but the evolution of jets with different param- 
eters will be similar to each other, if times are scaled to 
the time at which the gas begins to cool, and distances are 
scaled to the length of the jet at this time. For the cooling 
law, equation (2), the lengths will scale as 



Ic 



6 X lO'-" cm 



(l 



no 



> 1.5 cm~3 
and the times as 



2 X 1039 erg s-i 

-(6-2a)/(9-4c<) 



{3-2q)/{9-4q) 



(13) 



^cool 



■( 



no 



2 X 1039 erg s-i 

(7-2a)/(9-4a) 



(2-2a)/(9-4a) 



1.5 cm 



(14) 



If the jet power is doubled, the timescales will increase by 
22 per cent and the lengthscales by 30 per cent, while if the 
ambient density is doubled the lengthscales will decrease by 
35 percent and the timescales by nearly 40 per cent. 

Once the cooling becomes rapid at the head of the jet, 
when Ic — 7?.'^^'^icooi, the global evolution of the jet will en- 
ter a second self-similar regime with its evolution varying 
according to equation (^|), only with a rather smaller con- 
stant of proportionality. 



3.3 The hydrodynamic code 

We used the adaptive grid hydrodynamic code described by 
Biro et al. (1995) in axisymmetric mode. It solves the equa- 
tions of mass, momentum and energy conservation using a 
flux-vector-splitting scheme. The computation was carried 
out in a 5-level, binary adaptive grid. In our simulations we 
assume full ionization of the jet and ambient medium due to 
photoionization from the central UV source. We use the non- 
equilibrium cooling function described by Biro et al. (1995) 
with an additional term taking into account bremsstrahlung 
losses, Li,, given by 



Lb = 2.2<dxlQ~'^'' n^ripTi ergs^^ cm"^ 



(15) 



where n^ and Up are the electron and proton densities in 
particles per cm^, respectively, and Te is the electron tem- 
perature in Kelvin (Cox & Tucker 1969). We do not allow 
the gas temperature to drop below 10** K, where we assume 
it is maintained by the thermostat effect of [O ill] 5007-A 
emission. 

The computational domain was set to be 513x1025 com- 
putational cells and 5x1.25 times lO^^cm. The grid cells are 
smaller in radial direction than in axial direction for better 
resolution of the small distance scales of the radially ex- 
panding envelope. The initial jet radius was covered by 27 
cells in the axisymmetric simulation. The cooling distance 



5 was marginally resolved for Vs = 150 km s~^. The thick- 
ness dc of the cold envelope of swept up gas was near the 
resolution limit and was smeared out artificially over 2-3 
cells at shock speeds higher than 150 km s~^. This artifi- 
cially limits the compression with effects similar to those 
expected if magnetic fields in the ISM were present (with 
a magnetic parameter of the order of 5, although magnetic 
fields were not taken into account explicitly). However, in- 
creasing the resolution further would lead to prohibitively 
large computing times for these initial studies. Quantitative 
results will therefore be only order of magnitude estimates, 
but the qualitative results, as discussed in this paper, will 
not be significantly altered. The boundary conditions are 
reflective on the axis and on the left side of the computa- 
tional domain (except for the inflow condition where the jet 
is injected, see Fig. hi). The top and right boundaries have 
outflow conditions. 



3.4 Emission maps 

3.4-1 Optical emission 

We calculate the Hq emissivity e using radiative recombi- 
nation (Case B) following AUer (1984). In the temperature 
regime of our simulations (T > 10*K) only recombination 
contributes considerably to the emissivity. We assume full 
ionization of the hydrogen, which would mainly be due to 
photoionization from the UV source at the centre of the 
galaxy and, to some extend, from shock ionization in the 
bowshock region. 

e = 4.16 X 10-^^rgcm^s-i(Tl''«««10" °*2*/^^)-^ (16) 
where T4 is the temperature in units of lO'^K. 



3.4-2 Radio emission 

Only relativistic electrons (and perhaps positrons) con- 
tribute significantly to the radio emission of jets. Since it 
is still unclear how these high energy electrons are gener- 
ated (most probably through re-acceleration in shocks), we 
assume that this population of fast electrons shares the same 
hydrodynamics as the non-relativistic fluid of electrons and 
protons. A more detailed description would require the con- 
sideration of two jet populations to account for the non- 
relativistic gas (electrons and protons) with ratio of specific 
heats of 5/3, and the ultrarelativistic electrons with ratio 
of 4/3. Duncan, Hughes, Opperman (1996) have performed 
preliminary simulations of highly relativisitic jets using a 
variable adiabatic index and found that the rest frame vari- 
ables vary smoother along the jet axis, with smaller ampli- 
tude maxima at the bowshock, and in the vecinity of the 
contact surface. We ignore these effects in the present sim- 
ulations. 

In order to calculate the synchrotron emission from the 
jet whose hydrodynamics is modeled as above, we need to 
establish how the internal energy is distributed among the 
relativistic electrons. Observations suggest that this distri- 
bution follows the usual power law N{E) dE = NoE^^AE, 
with spectral index p, and energies above a minimum value 
Erain- Assumiug that the relativistic electron energy density 
U, and number density N are proportional to their corre- 
sponding thermal quantities calculated by the hydrodynam- 
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ical code as a function of position in the jet, the previous 
power law is determined by the equations (Gomez et al. 
1995) 



iV„ = [W(p-2)rMAA(p-l)] 



lP-2 



and 



En 



N p-1' 



(17) 



(18) 



The synchrotron emission is also determined by the 
magnetic field, and since we are neglecting its influence on 
the fluid dynamics, we assume that the magnetic energy den- 
sity remains a fixed fraction of the particle energy density, 
which leads to a field of magnitude proportional to W^'^^ . We 
refer the reader to Gomez et al. (1993, 1995) for a detailed 
discussion of the radio emission calculations. 



4 RESULTS AND DISCUSSION 
4.1 Structure 

We show the results of our simulation at a time when the 
jet has reached a distance of approximately 150 pc from the 
injection point. This is 127 000 years after the jet switched 
on. 

Fig. |l| shows the density (top) and the pressure (bottom) 
distribution on a logarithmic scale. We can identify several 
characteristic features which are marked in the schematic 
view shown in Fig. ^. The jet itself shows several recoUima- 
tion shocks on the axis, best seen in the pressure image. 
At the tip of the jet, the termination shock generates high 
pressure and high temperature jet plasma. This plasma flows 
back into the hot jet cocoon, which surrounds the high speed 
jet. The jet cocoon is surrounded by a thick layer of hot gas 
from the interstellar medium which has passed through the 
bowshock. Near the head of the jet we find a thin layer of hot 
shocked gas from the ISM which further back turns into an 
even thinner, but cold and very dense envelope of the jet co- 
coon. This envelope is what distinguishes the NLR jet from 
adiabatic high power extragalactic jets. The fact that this 
envelope starts at a significant distance from the front of the 
bowshock distinguishes NLR jets from the regime of strongly 
cooling galactic jets (e.g. Blondin, Fryxell and Konigl 1990). 

We find that the radius of the cold envelope in the simu- 
lation is approximately 7.8x lO^^cm and the expansion speed 
of the cocoon is somewhat above 100 km s~^. The peak den- 
sity in the envelope varies from values between 30 to 40 cm~'^ 
over most of its extent, up to 100 cm~"^ in the dense clouds. 
This corresponds to compression factor of 30 and 70, respec- 
tively. This is factors of 2 too low compared to Eq. (^) , which 
is probably due to numerical smearing. This efi'ect mimics 
a magnetic parameter of the order of 5 (see Eq.|l^), which 
corresponds to a ambient magnetic field near b^G. The in- 
tegrated Hq luminosity at this timestep is 5.8x lO^'^erg s~^ 
and the envelope has a thickness of dc ~ 1 — 2x10^* cm 
and volume of 14 ~ 2xl0^^cm'^. These values are in good 
agreement with our analytical estimates and the results from 
optical observations of Seyfert galaxies with associated NLR 
radio jets. 

We find that the non-stationary dynamics of jet propa- 
gation through interstellar medium have a strong infiuence 



on the structure of the cold shell. The quasi-periodic recoUi- 
mations of the jet cause the advance speed of the jet head to 
oscillate, thereby accumulating jet material in a restricted 
nose region for some period of time. It is then released into 
the cocoon almost explosively, when the jet expands and 
slows down again. This imposes a characteristic arc struc- 
ture onto the bowshock and the cold shell of swept up gas. 
These arcs cool almost as a unity once catastrophic cooling 
sets in. This can be seen in Fig.^, where the density is shown 
at three different times with the same separation between 
them. While the head of the jet has advanced considerably 
between the first and second panel, the point where the cold 
dense envelope starts has not advanced appreciably. How- 
ever, later, after the same timestep, the full arc has cooled. 
Therefore, a cooling distance for the bowshock measured 
from the head of the jet, as discussed by Taylor et al. (1989, 
1992), can at best be only an average property. 

At the positions where neighbouring arcs meet, the two 
shock waves cross at an angle near 90° . This interaction will 
produce a Mach shock at the vertex with an effective speed 
of \/2 times the speeds of the incident shocks (supposing 
that they are the same). Correspondingly, the compression 
and the emissivity of the gas increase by factors of ~ 2 
and ~ 4, respectively. As a result of the radiative nature of 
the crossing shocks they merge to form short high density 
filaments and clouds, marked as "high density spots" and 
"intruding dense filaments" in Fig.^ In the simulation the 
spots have peak densities between 80 and 100 cm~^, consis- 
tent with the previously estimated additional compression 
by factors around 2. We suggest that these spots could be 
identified with the brightest discrete clouds in the NLR of 
some Seyfert galaxies. Note that the spacing between these 
spots is directly related to the time dependent quasi-periodic 
recoUimation of the jet and therefore contains information 
about the history of the interaction between the jet and its 
environment. Several Seyfert galaxies show a knotty NLR 
structure and the distances between the strongest emission 
line clouds is indeed similar but rather smaller than the dis- 
tances between the radio knots (e.g. MKN 3, Capetti et 
al. 1995a; MKN 6, Capetti et al. 1995b; NGC 4151, Boksen- 
berg et al. 1995). 

At some distance from the jet head the arc structure in 
the envelope is lost and it becomes roughly cylindrical justi- 
fying the approximations of Section |^ (Fig.^. However, the 
high density clouds and filaments created at the intersection 
points of these arcs are retained. When the envelope frag- 
ments (as can be expected in a full three dimensional sim- 
ulation) these regions will be seen as discrete NLR clouds 
and filaments. Since the creation of these arcs and clouds 
is related to the recoUimations of the jet, we have found a 
direct relation between the optical structure of the NLR and 
the internal non-stationary shock structure of the jet. 

The expansion of the overpressured jet cocoon into the 
undisturbed interstellar medium is similar to a supernova 
explosion. Some theoretical aspects can therefore be treated 
in an analogous fashion. This applies in particular to the sta- 
bility of the shock and the postshock layer of swept up gas. 
The stability of radiative shock wave has been considered 
by a number authors with varying emphasis (e.g. Vishniac 
1983, Bertschinger 1986). The radiative shock instabilities 
discussed by Bertschinger (1986) are present in our simula- 
tions and produce the pressure variations near the cold en- 
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velope. These are likely to be disruptive and thus contribute 
to the formation of NLR clouds and filaments. From the 
results presented by Jun, Norman and Stone (1995) it can 
be expected that the shock region in our simulation is also 
Rayleigh- Taylor (R-T) unstable, with fastest growth of fea- 
tures on the scale of the thickness of the cold layer of swept 
up gas. The onset of these instabilities can be observed as 
pressure variations and small ripples in some sections of the 
cold envelope (Fig.^). As these instabilities grow, they might 
contribute to the fragmentation process of the envelope and 
form further NLR clouds, in addition to the high density 
clouds produce at the intersection of the arcs discussed in 
the previous section. However, magnetic fields have a stabi- 
lizing effect and therefore different combinations of ambient 
magnetic parameters and jet powers might lead to two cat- 
egories of evolved cocoon envelopes: R-T stable narrow line 
regions with fewer small clouds and R-T unstable ones with 
a high number of small NLR clouds with separations of a 
few times the thickness of the cold layer (i.e. ~ Ipc). More 
detailed investigations of this process will have to be carried 
out to quantify this prediction over a range of parameters in 
order to compare it with recent and future high resolution 
observations of the NLR in Seyfert galaxies. 

The Stromgren column given by Eq.(^ is similar to the 
thickness of the cocoon envelope and the size of the clouds. 
It may well be that the highest density clouds can develope 
neutral cloud cores. As a result, NLR clouds could be sites 
of jet induced star formation (Van Breugel and Dey 1993, 
Tresch-Fienberg et al. 1987). 



4.2 Emission 

In Fig.^ the Ha emissivity distribution is displayed on a 
logarithmic along with the radio emission superimposed as 
contours in the lower panel. As expected, most of the optical 
line emission comes from the thin envelope of cold postshock 
material. The Ha luminosity obtained from this axisymmet- 
ric simulation is 5.8 x lO^^ergs"^. This is well within our 
expectations considering the simplicity of our analytical es- 
timates. The length of the NLR jet in MKN 1066 is similar 
to that in our simulation, although the transverse extent 
could be 2-3 times smaller. From observations by Bower et 
al. (1995) we determined the Ha flux coming from the area 
covered by the NE jet in MKN 1066 to be 3.3 x 10^^ergs"\ 
which is in excellent agreement with our simulation, taking 
into account the the simplifications of the model and the 
smaller radius of the cocoon in MKN 1066. 

/_From our model it can be expected that higher power 
radio jets yield faster cocoon shocks, producing broader 
emission lines. This will result in a correlation between 
line width and radio power, which is consistent with the 
observed correlation between these quantities (Wilson and 
Willis 1980, Whittle 1985). Our initial study does, however, 
not allow a more quantitative conclusion about the correla- 
tion. An extensive survey of the parameter space and other 
factors which influence the line width, like galaxy rotation, 
is needed for a detailed theoretical study of this correlation 
in terms of our model. 

The calculated line spectrum is shown in Fig.^ for three 
different viewing angles. In the spectrum, emission from 
the background has not been included in order to focus on 
the expanding envelope. Inclusion of the background would 



produce a narrow peak at zero velocity. Within the vol- 
ume of this simulation the background emission amounts 
to 4.8 X lO^^ergs"^, an order of magnitude less than con- 
tained in the envelope. The strength of this central peak, 
compared to the broader contribution from the envelope, is 
a measure of the ratio of disturbed and undisturbed gas. 
Because of the cylindrical shape of the expanding envelope 
and the small velocity of the gas along the shell, the line 
shape of the envelope is double peaked for viewing angles 
far from the axis but becomes almost rectangular if the line 
of sight is near the jet axis. The line width of the expand- 
ing cylindrical shell is, of course, not only a function of the 
expansion speed (which is approximately the same as the 
shock speed), but depends also on the angle between the 
symmetry axis and the line of sight. This is in apparent dis- 
agreement with results from Whittle (1985), which suggest 
that there is no strong correlation between the orientation 
of the galaxy and the NLR line width, assuming that the 
jet follows the rotation axis of the galaxy. However, obser- 
vations of radiation cones (e.g. Wilson and Tsvetanov 1994), 
which at least can be expected to have a closer link to the 
radio jets, and point symmetric bending of jets in Seyferts 
(presumably as a result of the interaction with the rotat- 
ing environment, Wilson and Ulvestad 1982, Steffen et al. 
1996) suggest that the jets often are not aligned with the 
rotational axis of the galaxies. This would strongly weaken 
any correlation between line width and galaxy rotation. 

The contours in Figj^ show the radio emission at 10 GHz 
as calculated from our model described in Section 3.4.2| . A 
spectral index in the power law energy distribution of the 
electrons of 2.4, and magnetic field strength at the jet inlet 
of lOO/^G have been used. A total flux for the jet of 4 mjy 
is obtained. As expected, most of the radio emission comes 
from the hotspot region near the head of the jet, with some 
low brightness emission from the recoUimation shocks. There 
are a number of Seyfert jets in which the the head of the jet 
is dominant at radio wavelength. However, in most cases 
there are a number of strong radio knots along the jet. More 
so in the case of NGC 1066, which has rather continuous 
emission along the jet (at the resolution of the VLA) and no 
hot spot near the head. This highlights the problem of the 
nature of the knots seen along radio jets in general. Internal 
shocks can produce this kind of knots. In our simulations 
we find two different types of internal shocks: the recoUi- 
mation shocks and the shocks produced by the interaction 
with the turbulent cocoon near the head of the jet. From our 
simulations and our model for the synchrotron emission we 
find that these are not strong enough to explain the radio 
emission knots seen in most of the Seyfert jets. One possi- 
ble solution are time variations of the ejection parameters 
at the base of the jet. Gomez et al (1996), Hughes, Duncan, 
and Mioduszewski (1996) and Komissarov and Falle (1996) 
have shown that variations of the ejection velocity can pro- 
duce strong internal working surfaces which visible as strong 
radio knots in the jet. Such variations in the jet properties 
possibly have some infiuence on the structure of the cold co- 
coon envelope and thus on the NLR. In a forthcoming paper 
we will investigate the effects of variable jet properties on 
the NLR in more detail (Steffen et al, in preparation). 
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5 CONCLUSIONS 

Using hydrodynamic simulations with radiative cooling we 
confirm the basic structure of the interaction between a jet 
and the NLR in Seyfert galaxies as obtained in a geometric 
model by Taylor et al. (1989). However, we find new features 
of the interaction, like the formation of NLR clouds and a 
larger extent of the cold cocoon envelope. Due to the time- 
dependent propagation speed of the jet into the ISM, the 
cooling distance of the bowshock material as measured from 
the head the jet varies episodically with time. Our simu- 
lations show that there might be a link between the time- 
dependent internal recoUimation structure of the radio jets 
and the narrow line clouds. We expect that the distances 
between particularly bright knots along the symmetry axis 
are similar to the separation between the radio knots be- 
hind the main hotspot (assuming that the latter are due to 
internal recoUimation shocks). Estimates of the Stromgren 
column show that the dense clouds produced in this simula- 
tion could be partially neutral. They might represent cores 
for jet induced star formation in the NLR of Seyfert galaxies. 
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Figure 1. Density (top panel) and ijrcssurc (bottom panel) distribution on a logarithmic scale after 127 kyr. Note the high-density 
cocoon envelope after the swept up ISM-material has cooled to the equilibrium temperature of 10^ K. The recoUimation shocks and the 
very high pressure near the jet shock can clearly be seen in the pressure image. 
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Figure 2. A schematic summary of the features seen in our hydrodynamic simulations. It includes the jet with rccoUimation shocks, the 
jet-shock, the cocoon of jet plasma, the hot bowshock and the cold layer shocked ISM. 
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Figure 3. The bowshock structure (using the Ha emissivity) is shown at three different times with the same separation between them. 
While the head of the jet has advanced considerably between the first and second panel, the point where the cold envelope starts has 
not advanced appreciably. After a further timestep, a full arc has collapsed and cooled and has become very bright. 
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Figure 4. The Hq emissivity map is shown on a logarithmic scale along with the calculated radio emission at 10 GHz superimposed 
as contours (radio contours are drawn for 0.6, 5 and 50 percent of peak brightness). Most of the optical emission comes from the thin 
envelope of swept up and cooled ISM gas. Particularly high emissivity regions are found where the bowshock 'arcs' intersected when still 
near the head of the jet. The radio emission is concentrated near the head of the jet, where temperature and pressure of the jet plasma 
is highest after passing through the jet shock. 
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Figure 5. The spectral line shape of the Ha emissivity distribution as seen from different angles between the line of sight and the jet 
axis is shown. For large angles the line is double peaked, but for angles near the jet axis the shape is rather rectangular. 
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